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We report time-resolved x-ray scattering measurements of the transient structural response of 

■^1- , the sliding Qi charge-density-wave (CDW) in NbSe3 to a reversal of the driving electric field. The 

^\ ■ observed time scale characterizing this response at 70K varies from ~ 15 msec for driving fields near 

threshold to ~ 2 msec for fields well above threshold. The position and time-dependent strain of 

the CDW is analyzed in terms of a phenomenological equation of motion for the phase of the CDW 

order parameter. The value of the damping constant, 7 = (3.2 ± 0.7) x 10 -19 eV ■ seconds • A -3 , 

is in excellent agreement with the value determined from transport measurements. As the driving 

field approaches threshold from above, the line shape becomes bimodal, suggesting that the CDW 

Cn ' does not depin throughout the entire sample at one well-defined voltage. 
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The fundamental statistical physics describing systems which are so far out of equilibrium that the notion of a 
partition function is not valid is currently not well understood. Our goal is to study the structural response of 
very simple systems as they are driven between two distinct steady-state configurations, using time-resolved x-ray 
scattering. The particular experimental systems we chose to study are the charge-density-waves (CDWs) found 
in quasi-one-dimensional metals. The structure oLthe pinned CDW state has been measured at high resolution. 
Theory and experiment are in excellent agreementoa. On the other hand, the dynamics of CDW systems are less 
well understood, particularly the structural aspects. Although a large number of electronic transport experiments 
have been performed, the results of these experiments are difficult to interpret in terms at-microscopic models and 
time-resolved structural data are limited and restricted primarily to the K0.3MO3 systemrra 

In order to acquire the data, we developed a time-resolved, high-resolution x-ray scattering system capable of 
measuring the time-evolution of the structure of the sliding Qi charge-density- wave in NbSe3 as the direction of the 
driving electric field is reversed. We are able to interpret our data using an equation of motion for the phase of the 
CDW order parameter which pertains at low temperatures and at large applied fields. 

In the next section of this paper, we review the elementary physics describing CDWs in one-dimensional metals 
and discuss the structure of both stationary and sliding CDWs. In the third section, we review the standard theory 
describing x-ray scattering from CDW systems and develop some of the extensions necessary to describe evolving 
systems. We explicitly solve an equation of motion describing the wave number of the CDW. In the fourth section, we 
review the details of our NbSe3 samples and our experimental apparatus. In section five, we present our time-resolved 
x-ray scattering data. We conclude with a brief discussion of our results and an appendix of mathematical details. 

II. CHARGE-DENSITY WAVES 

In his 1955 solid state text, Peierls demonstrated that an ideal one-dimensional metal crystal is unstable to the 
formation of a CDW state at low temperatures^. In the CDW state, the conduction electron density, p c (x), is given 
by 

p(x) = p + p c cos(Qx + (j)) (1) 

where p is the mean conduction electron density, Q — Ikp is the CDW wave number, p c is the CDW condensate 
density and (j) is the phase. In the presence of the lattice-distortion wave, the position of the j th atom, Xj, is given by 



Xj = ja + u sin(Qja 



(2) 



where u is the amplitude of the distortion wave, and a is the real space lattice constant of the undistorted lattice. The 
conduction-electron-density wave and the lattice-distortion wave have the same phase, 4>. In the literature, CDWs 
are frequently characterized by a complex order parameter, A.e l ^ x K The amplitude, A, is usually taken to be the 
energy gap, which is linearly proportional to the amplitude of the lattice-distortion wave. The lattice-distortion wave 
and the associated conduction-electron-density wave are collectively referred to as a charge- density-wave. 

The properties of CDW systems have been studied extensively during the past 15 years. A large fraction-,of this 
research has been devoted to studies of the nonlinear electronic transport exhibited by some CDW systemsQ. This 
nonlinear transport originates in an incommensurate CDWs ability to slide rigidly through an ideal crystal lattice 
without frictionO. Real materials always contain lattice defects or impurities. These localized defects have energetic 
preferences for the phase and tend to "pin" the CDW in small regions, hence the term "defect pinning." 

Two phenomena are associated with defect pinning. First, the system can take advantage of the pinning energy 
of the defects by elastically distorting the CDW. This elastic distortion can be described by allowing the phase, 4>, 
to become a function of position. Clearly, this position dependence destroys the long-range periodicity of the CDW 
state. The phase-phase correlation function, < e l ^ Xl ) e~'"^ X2 ) > ~ e -l x i- x 2|/^ j s frequently used to describe this loss 
of long range periodicity and is characterized by the length scale £. A large fraction of the theory of and experiments 
performed on CDW systems have been devoted to studying this loss of long-range order. Second, and of particular 
relevance to this experiment, arbitrarily small fields are no longer rftapable of causing the CDW to slide. The CDW 
remains pinned until the applied field exceeds a threshold value, Ej€rcl. Experimentally, the onset of sliding is signaled 
by the presence of a non-ohmic contribution to the conductivity when the current exceeds the threshold current, It- 



The Sliding State, Boundary Conditions, and Phase Slip 



In any experimental realization, electrical contacts are required to produce the electric field which drives the CDW 
into the sliding state. The boundary conditions imposed by these electrical contacts play a crucial role in determining 
the structure of the sliding CDW state. To achieve a steady-state collective current due to the CDW sliding, current 
must be injected at one electrode and extracted at the other. Equivalently, CDW phase fronts must be added at one 
electrode and removed at the other. In Fig. [l|, the dashed lines indicate lines of constant phase of the CDW. In the 
region of the sample between the contacts, the CDW phase fronts are moving when the CDW is sliding. In the regions 
of the sample which are outside the electrical contacts, there is no collective current and the CDW phase fronts are 
stationary. 

This intuitive picture suggests a simple model for the current conversion process, originally proposed by J.C. Gill to 
explain the transients observed in his pulsed electronic transport measurementalj. Consider the contact on the top 
of Fig. pi. Before the electric field is applied, the CDW is uniform. As the electric field is applied, phase fronts in 
between the electrodes try to move down, dilating the CDW near the electrode on the top. This strain can be relieved 
by nucleating a dislocation loop, which quickly grows to be the size of the sample, inserting a new phase front as it 
does so. This mechanism of nucleating new lines of phase is called phase- slifiJJL. 

In electronic transport experiments, the phase-slip mechanism is studiedtJu in a "transposed" configuration. In 
the "normal" four-probe configuration, current is injected at-the .eiids of the sample and the voltage is measured at 
contacts which are separated by a distance L. MeasurementaiJ'tsliil find that, for a given Icdw, there is a difference 
between the voltage measured in the normal and transposed configurations. This difference is characterized by the 
phase-slip voltage, V ps . 

e-dimensional theory of the steady state current conversion process has been developed by Ramakrishna et 



ikS. 



In this theory, the dislocation loops are thermally nucleated in the presence of a background strain field. As 
in any nucleation and growth problem, loops above a critical size grow and those below the critical size collapseO. 
The resulting equations of motion for <f>(z, t) are linearized and solved for the time-independent component, <fii{z). 
The steady-state solution to this model has a parabolic profile and is zero at the electrodes, 
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Here, K z is the CDW elastic constant, and z is the distance from the midpoint between the contacts. The positional 
dependence of <f>\ creates a positional dependence, of Q, the CDW wave number. Using synchrotron-based high reso- 
lution x-ray scattering techniques, DiCarlo et a/O have observed a position, temperature and electric field dependent 
strain of the sliding Qi CDW in NbSe3 consistent with the specific predictions of this model. The goal of the present 
experiments is to measure the time dependence of this strain field as the sign of the electric field is reversed. 



III. X-RAY SCATTERING FROM CDWS 

X-ray scattering is a nearly ideal probe of the periodic structure of a CDW. In this section, we will review the 
standard theory relating the observed x-ray scattering to the microscopic structure of the CDW. Then, we will discuss 
some extensions necessary to describe time-dependent systems. 

Kinematic Scattering Theory 

In the kinematic approximation, the scattered intensity measured in an x-ray scattering experiment is [directly 
proportional to the spatial Fourier transform of the equal-time electronic density-density correlation functions 

S(q, t)ocf d*! J d X2 e *H*i-*»> (p(x 1 )p(x 2 )) p(t) . (4) 

In this expression, q is the scattering vector and p is the electronic charge density. The angular brackets indicate an 
ensemble average over all accessible configurations of the system. The subscript P(t) identifies the time-dependent 
probability distribution to be used when taking this average. 

At low temperatures, the electronic charge density of an ideal, monatomic crystal can be approximated by 

p(x) = ^p J 4(x-R j ) (5) 
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where pa(x) is the electronic charge density of an isolated atom and {Rj} is the set of equilibrium lattice sites. The 
sinusoidal lattice distortion associated with the CDW state can be introduced by letting 

Rj -*■ R, + u sin(Q • R^ + 0(Rj)) (6) 

We then proceed to calculate the T = structure factor, <S(q, t), arriving at the expressionEffiil 

5(q) ~ ^(^pJlJolq-^p^e^^ + lJxlq-u)! 2 ^^^-^^^^ )-^^)^.-.! (7) 

Here, |F(q)| 2 is the atomic form factor and J n {x) is the n th order Bessel function and the notation specifying P(t) has 
been suppressed. The first sum inside the square brackets produces the Bragg peaks associated with the undistorted 
crystal lattice. The second sum creates satellite peaks separated by ±Q from each Bragg peak. Since this is still 
a T = model, the angular brackets in this term indicate an average over different realizations of the quenched 
defect distribution. The line shape of the CDW satellite peaks is given by the spatial Fourier transform of the 
equal-time-phase-phase correlation function. 

In these measurements, q • u << 1 so we can use small argument expansions for the Bessel functions to simplify 
Eqn. 0, producing 

5(q) ~ |F(q)| 2 J^e^^ + (q-u) 2 ^e i(q±Q >^(e* (0 )-^ R ^) + ---l (8) 

which exhibits both the necessary Bragg peaks and the q 2 behavior commonly associated with disorder scattering at 
the CDW satellite positions. 

Note that the x-ray measurements are dominated by the localized atomic core electrons. These x-ray scattering 
experiments on CDW systems probe the structure of the lattice-distortion-wave component of the CDW state. The 
conduction-electron-density wave plays only an incidental role in determining the form of the observed scattering. 

X-ray Scattering Theory for Evolving Systems 

In principle, one would like to analyze the data using an expression containing the full <S(q, t). Note that, as defined 
in Eqn. 0, the time dependence of <S(q, t) originates in the time-dependence of the probability distribution P(t). Thus, 
in a time-dependent measurement, one is studying the evolution of P(t). Frequently, the Fokker-Planck equation is 



used for such studies. Such an analysis is not out of the question-ijor this system. For example, the elegant formalism 
developed to study the kinetics of first order phase transitionsE3 might be applied to this problem. However, the 
signal to noise ratio of our data is not yet sufficiently high to warrant such an analysis. Instead, as in the experiments 
of DiCarlo et alx3, wc will only seek to understand the physics describing the evolution of the strain of the CDW as 
a function of time and position on the sample. In essence, we will be studying the evolution of the first moment of 
the probability distribution, n ot the , full distribution function. 

As in the steady state case,c3E3ll3 we begin with a one-dimensional Langevin equation of motion for <fi(x, t) 

7^ = *.A»*|M + %LE{t)+ W0 + **.*)■ 0) 

Here 7 is a damping constant, Fi mp (<p) is the quenched random pinning force and r)(z, t) is the Langevin noise field. 
To study the kinetic problem, we now follow a similar line of reasoning to that used by Ramakrishna et alS3 for the 
steady-state case. We begin by defining (j)(z, t) — 4>o(z, t) -\- <f>\(z) . Where, <fo represents the bulk solution to the sliding 
CDW problem and <j>\ is a small time-independent correction which describes the background strain field required 
to drive dislocation loop nucleation at the current contacts. For the transient problem, instead of requiring (f>i to 
be constant in time, we allow <f>i to have a slow time dependence. By "slow" we mean that <p$ fluctuates so rapidly 
relative to <p\ that 4>o is always in equilibrium with the current value of (j)\. We show below that the characteristic 
time scale for changes in <p\ is on the order of milliseconds, much-slower than any electronic time scales in this system. 
We then continue to use the same logic as Ramakrishna et aZ.E3il3, to develop an equation of motion for (f>i(z, t). In 
the limit where the driving electric field is far enough above threshold that the pinning force and the noise term are 
both negligible, this equation of motion becomes a simple driven diffusion equation, 

^ - *•*£ + (£) ¥>■ « 

We are interested in the solution to Eq. [10|for the simple case where the sign of the driving field, V ps (t), is reversed 
at t = and </>i(±-|) = 0. Laplace transform techniques produce a series solution which converges rapidly at early 
times. Conversely, the separation of variables technique produces a series solution which converges rapidly at late 
times. The details of the latter series solution are given in the appendix. Using the appropriate solution at a given 
time, t, we have a convenient representation of the solution for all times t > 0. 

IV. EXPERIMENTAL 

NbSe 3 

We chose to use NbSe3 for our experiments for a number of reasons. First, although a large number materials are 
known to exhibit CDW states, the CDWs in NbSe3 both exhibit sliding mode conduction and have been extensively 
studied in electronic transport measurements. Second, we are able to grow crystallographically perfect, single crystal 
whiskers of NbSe3. These structural defect- free samples are essential if one wishes to measure these small strainsE3. 
Third, in NbSe3 the Fermi surface is not completely gapped below the first Peierls transition. The system remains 
metallic and the electronic transport is still dominated by the normal carriers. Consequently the electric field must 
remain uniform inside the whisker, eliminating spurious polarization effects due to poor electrical contacts in regions 
not too close to the electrodes. Finally, the technology required to produce good electrical contacts has been well 
developed for NbSe3. 

On a microscopic level, the crystal structure of NbSe3 is monoclinic, with lattice constants a = 10.009A, b — 3.4805A, 
c = 15.629A, and (3 — 109.47°E-i Macroscopically, these crystals grow as long whiskers with a rectangular cross section. 
The width of a whisker is typically ten times the thickness. Typical dimensions of a sample used in our experiments 
are on the order of 2 /im x 20 /im x 10 mm. Both the crystallographic b and b* directions are oriented along the 
whisker axis. The smallest macroscopic sample dimension usually corresponds to the a* direction. 

Two independent CDWs form in NbSe3. The Qi CDW forms at a temperature of Tp ± = 145 K and has the wave 
vector qi = (0 Q\ 0) where Q\ is is slightly temperature dependent and approximately equal to 0.243E3. The Q2 
CDW forms at roughly Tp 2 = 59 K and has the wave vector q2 = (0.5 Q2 0.5) where Q2 is approximately equal to 
0.263. All of the measurements reported in this paper were performed at temperatures above Tp 2 and below Tp 1 . 



Experimental Configuration 

The x-ray scattering measurements were performed at the F2 experimental station at the Cornell High Energy 
Synchrotron Source (CHESS). The storage ring was running at an energy of 5 GeV and the stored positron current 
typically decayed from 80 to 40 mA during a 50 minute fill cycle. For the purposes of comparing different data sets, 
we have normalized the data to counts per second at 100 mA ring current. A Si(lll) double-bounce monochromator 
selected a wavelength of 1.5A from the white x-ray beam produced by the 24-pole wiggler. A flat Au-coated mirror in 
the hutch suppressed harmonics of the fundamental wavelength passed by the monochromator. The sagittally-bent 
second monochromator crystal focused the x-ray beam in the out-of-scattering-plane direction at the sample position. 
Tantalum slits restricted the x-ray spot size at the sample to approximately 0.8 mm x 3 mm. The resulting x-ray 
beam contained 4 x 10 10 x-rays/second/100 mA of stored positron current. The scattered x-rays were analyzed by a 
triple-bounce channel cut Si(lll) crystal and detected by a standard Nal(Tl) scintillator and photomultiplier tube. 

The NbSe3 whiskers were mounted across a 4mm hole in an alumina substrate using silver paint, which also provided 
the electrical contacts. The distance between the electrical contacts was L = 5 mm. In these measurements the sample 
was oriented such that the length of the sample in the beam was 0.8 mm, centered at roughly the | position between 
the contacts. r— ■ 

To study the transient structural response of the CDW, we made a stroboscopic measurements!! The CDW is 
subjected to a continuous square wave voltage wave form. The half period of the square wave is divided into time 
intervals of equal size. X-rays detected during a particular time interval are summed over a large number of square 
wave periods. Typical count rates at the CDW satellite were on the order of 50 counts/second @ 100 mA of ring 
current. Therefore, in order to obtain reasonable counting statistics in a 4 msec wide time bin, we sum over roughly 
10,000 voltage wave form cycles. Transient structural effects due to ohmic heating of the sample were eliminated by 
switching only the direction of a constant magnitude current, thus keeping the ohmic heating constant. Time-resolved 
measurements of the (0 2 0) Bragg peak as the current direction is switched do not exhibit any measurable change, 
ruling out transients due to changes in the lattice (e.g., a strain of the lattice due to a bend of the whisker driven by 
the electric field). 

All the data presented in this paper were taken on the (0 1 + Qi 0) CDW satellite. The (0 ± (1 + Q\) 0) 
CDW satellites are particularly attractive for high resolution work because they are nearly non-dispersive for Si(lll) 
monochromator/analyzer optics and because the (0 1 0) structure factor is zero in NbSe3, eliminating any diffuse 
background from the 4j- tail of the Bragg peak. The measured q-space resolution (FWHM) of our diffractometer in 

this region of reciprocal space was Sq t = 7.6 x 10 _5 A _1 , Sqi = 2.7 x 10~ 4 A _1 , and Sq z ~ 1 X 10~ 2 A _1 . 

V. DATA 

All of the data shown in the figures were taken at 70 K. Using the standard technique, the depinning current 
was measured by using a lock-in amplifier to measure the differential resistance as a function of current. At 90 K, 
It = 10 ± 1 mA. Based on this measurement, we estimatecH that at 70 K, It — 20 ± 3 mA. 

Figure illustrates a typical data set: the observed intensity, I(q,t), is plotted as a two-dimensional surface. The 
same data set is also shown in Fig. 13 as a contour plot. The evolution of the CDW satellite as the system is driven 
by the voltage wave form shown in the inset is clearly visible in either figure. The solid lines in Fig. @ are contours 
of constant scattered intensity. The contours are separated by 30 counts/second @ 100 mA. The magnitude of the 
driving current wave form in this case was 21 mA, slightly greater than the threshold current, It — 20 mA. After the 
current direction is switched, a slow shift in the position of the CDW satellite is visible with a time scale on the order 
of a few milliseconds. The period of the driving wave form is chosen so that the CDW satellite peak stabilizes at its 
steady-state position before the next switch in the direction of the current. 

Close inspection of Figs. or reveals that the satellite peak does not have the same shape or width for the different 
current directions. The data exhibits several other asymmetric features. In particular, the position and width of the 
line shape of the zero-field cooled state differ from those of the pinned state which a sliding CDW relaxes into after 
the driving field is removed. The details of these asymmetries are also position dependent. The asymmetry of the 
line shape of the driven sliding state is more clearly seen at lower driving fields. For example, Fig. shows a series of 
longitudinal scans through the CDW satellite from different time intervals as the current direction is switched. These 
data were taken at 70 K and with a driving current of 17 mA. Each scan is offset by 35 cps @ 100 mA for clarity. The 
top and bottom scans represent the steady-state for the two current directions. The remaining scans show that the 
CDW satellite peak evolves from a single peak to a bimodal line shape. Data sets taken with higher driving current 
magnitudes exhibited the same systematic structure. In general, as the driving current magnitude was increased, the 



area under the stationary peak decreased. At 29 mA, the stationary peak is no longer measurable. The width of the 
moving peak is slightly broader in the transient regions. 

To characterize the data, we fit to a two component line shape. One peak was constrained to have a time-invariant 
position, width, and integrated intensity. The integrated intensity of the second peak was also held constant, but the 
position and width were allowed to vary with time. Empirically, each of the peaks is well described by a pseudo-Voigt 
line shapeRa. The solid lines shown in Fig. |4| represent fits to the sum of two pseudo-Voigt functions. The values of 
X 2 for the fits were between 1 and 4 where we have assumed that the errors are due to counting statistics. 

Figure |)| displays the fraction of the total integrated intensity under the stationary peak as a function of the 
magnitude of the driving current. The fraction of the sample which remains pinned decreases monotonically as the 
current is raised. At a current magnitude of 29 mA, this fraction drops to zero. Due to systematic changes in the 
scattering geometry as the diffractometer was raised and lowered in order to vary z, we cannot directly compare the 
intensities observed in different sets. 

Figure @ shows the difference in position of the sliding and pinned peaks as a function of time for one complete 
period of the driving current wave form. The particular data set shown was taken at a .current of 29 mA but is 
representative of the data taken at all currents. Error bars represent standard fitting errors!!! 

Analysis and Discussion 

A simple interpretation of this data is that at driving fields near the threshold to sliding the CDW breaks into 
domains. Although a collective current is flowing, some portions of the sample remain unstrained and hence are 
presumed to be pinned. One possible explanation of the physics is that near an impurity site, the applied electric field 
may be insufficient to cause the CDW to slide. The same singularity at 2kp in the dielectric function, e(q, u> = 0), 
which drives the CDW phase transition alsp-causes the system to respond to a point impurity by creating a localized 
electronic state with wavenumber Q = 2fcp EH Indeed, this second effect is the origin of impurity pinning. However, the 
pinned portions of the sample should not exhibit the time-dependent strain of the sliding portions. If the domains are 
small enough, the CDW satellite peak should broaden and we do observe significant broadening of the CDW satellite 
peak in a* direction near threshold. Since the peak moves in the b* direction, it was not possible to characterize this 
broadening quantitatively. The peak width in the b* direction is roughly 25% greater in the sliding state than in the 
zero field cooled state. The observation that raising the current reduces the area under the stationary peak provides 
further support for this picture since at large enough fields (currents), the whole sample should be in the sliding state. 
In this scenario, the depinning current measured using the lock-in amplifier cannot be interpreted simply as the single 
current at which the CDW is depinned over the entire crystal. 

Quantitative analysis of these time-dependent strain results is somewhat problematic. Specifically, Eqn. |l0| applies 
only at driving fields well above threshold. The highest current at which we were able to obtain data was on the 
order of 1.5 x It- The model specifically predicts that the time constants for different directions should be the same; 
yet, the data are clearly asymmetric. Furthermore, there are not enough time bins in the transition region to make a 
reliable test of Eqn. [To] and hence the model. 

On the other hand, these data are sufficient to demonstrate that x-ray measurements produce results which are 
consistent with electronic transport measurements. To characterize quantitatively the time scales of the observed 
transients, we fit the relative peak positions to 

Q{z, t) = Q ± SQ (2e- f / ro - l) sin (H) (11) 

where r = a „ 7 . 2 . This equation describing the transient structural response is the derivative of Eqn. EG in the 
appendix. The solid line in Fig. |6| is the best fit to Eq. |ll|. The amplitude, SQ, and the nominal zero field value, Qo, 
were forced to be consistent for the current switch in both directions but the time constant To was allowed to differ. 
For the I = 29 mA data set shown in Fig. y, the time constant To obtained for the switch to larger values of Q was 
2.6 ± 0.4 msec and for the switch to smaller values of Q was 4.7 ± 0.7 msec. Using the latter value of To, Eqn. |2^ 
and K Z A 2 = (1.7 ± 0.25) x 10~ 2 eV ■ A' 1 £3 the value of 7 becomes (3.2 ± 0.7) x 10~ 19 eV • seconds • A" 3 . 

This value of 7 can be compared directly to the value of 7 obtained from the high-field CDW conductance using 
the relation 

— V 1 (12) 

2k F J (Jcdw 

Using the values ucdw = 2 x 10 5 fi _1 m _1 and p c = 3.8 x 10 27 m -3 appropriate for NbSe3, we find that 7 = 5.7 x 10~ 19 
eV • s • A -3 , in excellent agreement with the x-ray estimate. 



The time constants obtained as a function of driving current amplitude are shown in Fig. [7[ Error bars again 
represent fitting errors. For current magnitudes > 17 mA there is a monotonic decrease in To as the current magnitude 
increases. For 1=14 mA, however, there is a jump to shorter time scales, comparable to those observed at 29 mA. 
Similarly, the magnitude of the strain is much smaller. We do not have an explanation for this behavior but observe 
that this jump occurs at roughly the depinning current. It is possible that at currents -below the depinning current, 
the CDW is "sloshing" back and forth slightly as hypothesized by Feinberg and FreideEj. 

VI. SUMMARY & CONCLUSIONS 

In summary, we have developed a time-resolved high-resolution x-ray scattering system capable of measuring the 
transient structural response of a simple quasi-one-dimensional CDW system as it is driven between two distinct 
steady state configurations. We find that a simple Langevin equation of motion for the phase of the CDW, which is 
based on the notion of a background strain which is driving the nucleation of dislocation loops, gives a value of the 
phenomenological damping constant, 7, for the Qi CDW in NbSe3 of (3.2 ± 0.7) x 1CP 19 eV • seconds • A~ 3 . This 
value of 7 is in excellent agreement with the value obtained from electronic transport experiments. As the driving 
field approaches threshold from above, the line shape becomes bimodal and the area under the stationary component 
of the line shape increases monotonically. This suggests that, the CDW did not depin throughout the entire sample 
at one well-defined voltage. Clearly much more work needs to be done. In particular, experiments addressing the 
asymmetry in the response of the sample and simultaneous electronic transport measurements need to be performed. 
However, these preliminary x-ray measurements demonstrate the potential of the technique for studying kinetics in 
CDW systems. 
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APPENDIX I 

A series solution to Eqn. QO^can be obtained by performing a solution similar to a separation of variables calculation. 
Equation [H] has the generic form 

<f> t = a 2 (j) zz + w(z,t). (13) 

We begin by assuming a series solution of the form 

Mz,t) = JTa n (t) cos ( {2n+ L 1)7TZ ) (14) 

71 = ^ ' 

which both spans the space of functions for — -j < z < 4 and satisfys the boundary conditions. This solution happens, 
in this special case, to have the form of a conventional Fourier series. 

Substituting this series solution into Eq. U3 leads to a differential equation for the a n (t), 

9 f(2n+ 1) 2 tt 2 \ , , , , 

a' n {t) = -a 2 l± jpi Ja n (t)+w n (t) (15) 

where w n (t) is the coefficient of the expansion of the driving function w(z, t); i.e., 

w(z, t) = JT w n (t) cos ( (2n + 1)7rZ ) . (16) 



L 

»=o x 



The solution of Eqn. na is 



a n (t) = cxp 



(2n + l) 2 7r 2 a 2 ' 



exp 



(2n+l) 2 7r 2 a 2 ■ 

-T 



w n (r)dT + c n > (17) 



L 2 
To determine the constants of integration c„, we now use the fact that <f>(z, 0) = f(z); i.e., 

/(*) = £ MO) cos (<?!_+_)__:) ( l 8) 

so that the c n = cv„(0) are simply the expansion coefficients of the initial state 4>{z,Q). For our problem, w(z, t) = 
( k A 2 ) ( ff^ ) ( "P ) an d the initial state is given by Eq. [| Therefore, 

1 \fe Pc \ 4V ps L(-l) n (ig) 



i^ 2 A 2 y \2k F J 7r 3 (2n+l) 3 
Thus, using these {c„}, the solution for cf>i(z,t) is 

(2n + l)nz 
c n I ze ' - — i cos I ■ 

.1=0 ^ 

(2e-*/ T » - l) cos (H) (21) 



l(M ) = 2 c™ [2e-V- - l] cos i-__+-J!_! (20 ) 



~ c 
where 



(22) 



(2ra + 1) 2 tt 2 K Z A 
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FIG. 1. Schematic illustration of the experimental configuration. Normal current is injected and extracted from the current 
electrodes. CDW or collective mode conduction takes place only in the region between the current electrodes. 

FIG. 2. Scattered intensity measured near the CDW satellite as a function of the scattering vector, [0 Q 0], and time, t, 
plotted as a two-dimensional surface. This data set was taken at a temperature of 70 K and current of 21 mA. It — 20 ± 3 mA. 
The signal waveform, shown in the inset to Fig. 3, switches polarity at and 80 msec. The time bins are each 8 msec wide. 

FIG. 3. A contour plot of the same data shown in Fig. 2. The contours are lines of constant scattered intensity. Neighboring 
contours differ by 30 cps @ 100mA. The highest contour represents a total accumulation of 300 counts in a time bin. The signal 
waveform, shown in the inset, switches polarity at and 80 msec. The time bins are each 8 msec wide. 

FIG. 4. A succession of (0 Q 0) reciprocal space scans taken from different time slices of a data set taken at 17 mA. Each 
time bin is 0.004 seconds wide. The centers of the time bins, measured from the rising edge of the driving wave form, are 0.014, 
0.030, 0.050, 0.062, and 0.078 seconds. Successive scans are displaced by 35 cps @ 100mA. The solid lines are the best fit to 
the sum of two pseudo-Voigt functions, one stationary and the other moving. 

FIG. 5. The ratio of the area under the stationary peak to the total area under the peaks is plotted as a function of current 
at 70 K. These areas are determined using the two component pseudo-Voigt line-shape demonstrated in Fig. 4. The error bars 
represent standard fitting errors. 

FIG. 6. Positions of the moving peak are shown for a data set taken at a current of 29 mA (It — 20 ± 3 mA). The solid 
lines are the best fit to Eq. 11 in the text. 

FIG. 7. The characteristic time constants, To, derived from fits such as those shown in Fig. 6 are shown as a function of 
current. Empirically, the time constants to approach the different current directions are not the same. For / > 20 mA, the 
values of tq decrease monotonically with increasing current. 
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